arXiv:hep-ph/0112209v3 5 Feb 2002 


A Way to Reopen the Window 
for Electroweak Baryogenesis 

G. Servant 

February 1, 2008 


Enrico Fermi Institute, University of Chicago, Chicago, IL 60637 
High Energy Physics Division, Argonne National Laboratory, Argonne, IL 60439. 


Abstract 

We reanalyse the sphaleron bound of electroweak baryogenesis when allowing devi¬ 
ations to the Friedmann equation. These modifications are well motivated in the 
context of brane cosmology where they appear without being in conflict with major 
experimental constraints on four-dimensional gravity. While suppressed at the time 
of nucleosynthesis, these corrections can dominate at the time of the electroweak 
phase transition and in certain cases provide the amount of expansion needed to 
freeze out the baryon asymmetry without requiring a strongly first order phase 
transition. The sphaleron bound is substantially weakened and can even disappear 
so that the constraints on the higgs and stop masses do not apply anymore. Such 
modification of cosmology at early times therefore reopens the parameter space al¬ 
lowing electroweak baryogenesis which had been reduced substantially given the 
new bound on the higgs mass imposed by LEP. In contrast with previous attempts 
to turn around the sphaleron bound using alternative cosmologies, we are still con¬ 
sidering that the electroweak phase transition takes place in a radiation dominated 
universe. The universe is expanding fast because of the modification of the Fried¬ 
mann equation itself without the need for a scalar field and therefore evading the 
problem of the decay of this scalar field after the completion of the phase transition 
and the risk that its release of entropy dilutes the baryon asymmetry produced at 
the transition. 
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1 Introduction 


Since the appearance of the article by Kuzmin, Rubakov and Shaposhnikov in 1985 [|l| who 
suggested that baryon number violation at high temperature in the Standard Model might 
have played a role in the early universe and set the basis for the theory of electroweak 
baryogenesis, there have been extensive studies which explored with great scrutiny the 
details of the mechanism and its viabiblity both in the Standard Model and in the Minimal 
Supersymmetric Standard Model. The fact that the creation of the baryon asymmetry 
of the universe is in principle possible in the standard electroweak theory is indeed very 
appealing since it relies only on physics accessible at present colliders in contrast with 
GUT or Affleck-Dine scenarios. Electroweak physics is not though the only ingredient. 
Electroweak baryogenesis makes use of an intricate interplay between particle physics 
and cosmology which is much more involved than in the theory of nucleosynthesis. Let us 
summarize very briefly the main idea: In a universe with conserved vanishing B — L charge 
and standard Friedmann-Robertson-Walker (FRW) cosmology, the baryon asymmetry of 
the universe (BAU) can only be produced at the electroweak phase transition (EWPT). To 
fulhll Sakharov’s out-of-equilibrium requirement for baryogenesis [0 the phase transition 
has to be hrst order. In that case, bubbles of non zero higgs held vev nucleate from the 
symmetric vacuum and as they expand, particles in the plasma interact with the moving 
phase interface (bubble wall) in a CP violating way producing a CP asymmetric hux 
of particles into the symmetric phase. This CP asymmetry is converted into a baryon 
asymmetry by baryon number violating sphaleron processes in front of the wall. As 
the phase interface propagates, the baryons enter into the broken phase where baryon 
violation is inefficient and therefore baryon asymmetry frozen out |^. 

To compute the resulting net baryon asymmetry produced during the electroweak 
phase transition one needs to know: (1) The details of the non trivial baryon number 
violating processes at high temperature (B violation is achieved thanks to sphaleron tran¬ 
sitions interpolating between two neighboring N-vacua of the electroweak gauge theory). 

(2) The amount of departure from equilibrium at the electroweak phase transition. This 
necessitates an accurate calculation of the higgs effective potential at high temperature. 

(3) The generation and propagation of CP-violating particle fluxes. (1) is now well un¬ 
der control. The rate of B-violating sphaleron processes at high temperature has been 
performed accurately and checked independently using analytic perturbative and non per¬ 
turbative techniques both in the broken phase |Q-§[ and in the symmetric phase []^,^. 
(2) has also been quite well understood. It is at this point that cosmology enters into 
the game when comparing B violation rate with the expansion rate of the universe. One 
of the strongest constraint in electroweak baryogenesis comes from the requirement that 
baryons produced at the bubble wall are not washed out by sphaleron processes after they 
enter the broken phase. Imposing that sphaleron processes are sufficiently suppressed 
in the broken phase at the critical temperature leads to the so-called sphaleron bound: 
{(j){Tc))/Tc > 1. (0(Tc)) is the order parameter of the phase transition, the vev of the 
higgs field at the critical temperature. The derivation of this bound uses the expression 
of the expansion rate of the universe H given by the Friedmann equation evaluated in a 
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radiation dominated universe. {(j){Tc))/Tc on the other hand is given by the analysis of the 
effective scalar potential at high temperature and depends sensitively on the higgs mass 
so that the sphaleron bound translates into a constraint on the higgs mass. It has been 
shown that in the Standard Model the phase transition is not hrst order but just a cross 
over for a higgs mass above ~ 72 GeV Q so that electroweak baryogenesis is excluded 
given the present experimental lower limit of the Standard Model higgs mass. Another 
reason in favour of that conclusion is related to our third point (3): CP violation in the 
Standard Model is not enough to account for the baryon asymmetry of the universe . 


Consequently, efforts have concentrated on the Minimal Supersymmetric Standard Model 
which provides not only the possibility of achieving a first order phase transition but 
also additional sources of CP violation |^. While there is a common agreement on the 
identihcation of the major CP violating sources-the complex parameters /i and M 2 in the 
mass matrix of the charginos- the result of the calculation of CP violating currents on 
the other hand remains an open question [1^-1^. We will not consider this issue here 


and assume that there are sufficient sources of CP violation in the MSSM to account for 
the baryon asymmetry of the universe so that the strongest constraint on electroweak 
baryogenesis essentially comes from the sphaleron bound. 

We will therefore focus now on point (2) and on the condition {0(Tc))/Tc > 1. 
(0(Tc))/Tc measures the strength of the phase transition and depends crucially on the 
presence of a light stop tji (with small left-right mass mixing parameter At). It is difficult 
to compute since it is radiatively induced and perturbation theory is expected to break 
down at high temperature. Nevertheless, the validity of perturbation theory and results 
at two loop level [|^|^ have been conhrmed by lattice calculation [|T^ which enhance 
(0(Tc))/Tc by 10 — 15% compared with perturbative calculation. From the last results 
one can infer the constraints imposed by the sphaleron bound: 


IS 


110 GeV <mh < 115 GeV and 105 GeV < 165 GeV (1) 

As the LEP bound on the higgs mass has been pushed up to rrih ^115 GeV, the window 
(2) for electroweak baryogenesis has been seriously reduced. A large value for the higgs 
mass requires in particular a large stop mixing parameter At. On the other hand a large 
At makes the phase transition weaker and leads to values of {(j){Tc))/Tc in conflict with 
the condition of preservation of the baryon asymmetry. 

In summary, given the last bound on the higgs mass, electroweak baryogenesis in the 
MSSM starts to be in bad condition. However, if it turns out that the window i) is 
ruled out by experiment shall we definitely claim the death of electroweak baryogenesis? 
Adding new particle physics ingredients to enhance the phase transition would start to be 
arbitrary and electroweak baryogenesis would lose its main advantage over other baryo¬ 
genesis mechanisms that it relies on a well motivated and well defined particle physics 
framework. If there is no reasonable way out on the particle physics side and if one still 
wants to believe that the baryon asymmetry of the universe was created at the electroweak 
phase transition the only way to turn around the sphaleron bound is to modify cosmology 
and question the common assumption made on the thermal history of the universe. There 
have been only a few attempts in the past to study the possibility of relaxing the sphaleron 
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bound using alternative cosmologies p0|-[22 
in 


An interesting discussion has been offered 
n[] . All|^ these works have in common the fact that they modify the usual assumption 


of a radiation dominated universe by assuming the existence of a scalar held driving the 
fast expansion of the universe at the EWPT. One drawback with this approach is that the 
decay of this scalar held after the completion of the phase transition may release entropy 
and dilute the baryon asymmetry produced at the transition. 

In this paper we take a diherent point of view and discuss a new cosmological solution 
to weaken the sphaleron bound. We still assume radiation domination at the EWPT. 
Instead, what we modify is the Friedmann equation itself in a way motivated by results 
in brane cosmology. While the usual expansion rate is recovered at late times in particu¬ 
lar at nucleosynthesis, the expansion can be accelerated at high temperature with major 

m 


consequences for electroweak baryogenesis as already suggested in 


In section II we 


present the derivation of the sphaleron bound and discuss how modihcations of common 
assumptions on the cosmological evolution could substantially weaken that bound. In 
section III we motivate our parametrization of deviations to Friedmann equation by re¬ 
viewing generic aspects of brane cosmology and discuss under which general conditions 
these deviations may play a significant role at the time of the EWPT. Section IV makes 
some comments on the cosmology with a non standard Friedmann equation. In section V 
we exhibit our numerical results describing the parameter space allowing the preservation 
of the baryon asymmetry. 


2 Derivation of the sphaleron bound 


To derive the sphaleron bound, one starts with the “master equation” of electroweak baryo¬ 
genesis which gives the rate at which baryon number relaxes to its equilibrium value in 
the presence of fermions in the hot plasma 


dt 


= -VB{t)n\ 


com 

B 


( 2 ) 


VB{t) = M 


T3 


(3) 


is the density of baryons per comoving volume in an expanding universe, = 
a{t)^nB, where a{t) is the scale factor of the universe. VB{t) is the rate of the baryon 
number non conserving processes. Af is related to the number of fermionic degrees of 
freedom at thermal equilibrium which bias the free energy of the SU{2) vacuum. In 
the Standard Model, A4 = 13Nf/2, if one assumes that all fermionic degrees of freedom 
(including cr) are at thermal equilibrium. In the MSSM, Ai = 15Nf/2 The 


sphaleron rate, P^, is the rate per unit time and unit volume of fluctuations with changing 


^For exhaustiveness, let us mention another type of alternative: If the EWPT was preceded by a 
color and charge breaking phase, the ratio (dfTcD/Tc could be much larger than in the usual scenario. 
However, it has been shown that such two-stage EWPT is not cosmologically viable ||^. 
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of the topological number. In the broken phase, a perturbative analysis gives P,p6H28 




where 


C(r) = E,p(T)/r and 77 = 2.8 x lO'" 


kB 


( 4 ) 


(5) 


EgpiT) = ^J~^B(j)(Tc) is the energy of the sphaleron conhguration in other words 

the energy barrier between two adjacent topologically non equivalent vacua of the SU(2) 
broken phase. B{X/g‘^) is evaluated numerically and is a monotonic slowly varying func¬ 
tion of the higgs mass: In the Standard Model, B{0) = 1.5 —>• B{oo) = 2.7 [p0|,|^ 
value of K has originally been assumed to lie in the range 10’ ^ ^ ^ 10 
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The 
how¬ 
ever, recent non perturbative calculations ^ have shown that it is reasonable to use the 
perturbative estimate (|[) provided that n rsj 0(1). 

Integrating (|^) between the time tc at which the electroweak phase transition takes 
place and some time t leads to 


r^COm( + \ ni 


( 6 ) 


where we have introduced the dilution factor S. The most efficient dilution of the baryon 
asymmetry takes place just after the phase transition at a temperature near the critical 
temperature Tc so that it is not a bad approximation (we will check it shortly) to replace 
the integral in (P) by the value of the integrand at T = Tc. 

where At is the time interval corresponding to a range of temperatures AT ~ — T 2 = 

(130 — 110) GeV, T 2 being the temperature at which dilution is inoperative. We note that 


n 


com 
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^3 _ com 

(X Tl — S 


Ub 


( 8 ) 


where s™™ = a^s is the entropy per comoving volume. In a radiation dominated universe 
s oc g^{T)T^. At this stage, we need the relation between time and temperature and that 
is where cosmology comes in. However, we want to remain as general as possible and will 
not assume any specihc form of the Friedmann equation for the moment. We will only 
assume that the universe is expanding adiabatically, d{pa?) = —pd{a?) {p is the pressure), 
and that the universe is radiation dominated, leading to the relationship between the 
energy density p and the scale factor of the universe a(t): p{t) = pQ^ait)/aQ)~^ where p is 
given by 


P 



(9) 
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5'* counts the number of degrees of freedom at thermal equilibrium. At T ~ 100 GeV, 
( 7 * = 106.75 in the Standard Model. From (|^) we get a{t) oc oc g7^''^{T)T~^ so 

that s™™ scales as oc gl^‘^{T). We can rewrite the dilution as 


?(t) ( a.(T) A ‘-'‘‘ 
“fe) \3.(Tc)) 


( 10 ) 


hb/s measures the baryon asymmetry. The value of the asymmetry today deduced from 
nucleosynthesis is 1.7 10“^^ <8.9 10“^^ so that 
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( 11 ) 


The value of the baryon asymmetry which can be produced at the EWPT, —(G), is still 
an open question mg , however, it is generically difficult to produce a large asymmetry 
and one can reasonably write the bound ns(G)/s(A) ^ 10“® which means that for the 
baryon asymmetry not to be washed out, one demands S > 10“®. This leads to 


In + 71nC(T,) - ((B) + ln(r,A() ,J0 (12) 

Note that equation (^) is general and does not depend on the precise form of the Fried¬ 
mann equation. It only relies on the adiabaticity of the expansion and on the radiation 
domination assumptions. It is when replacing At by its expression in terms of the tem¬ 
perature that we make use of a particular Friedmann equation. Let us now express this 
condition using the standard Friedmann equation: 


H = 



P(T) 



(13) 


mpi = Mpi/\/^ is the reduced Planck mass in four dimensions. This leads to 


p{t) = po 1 + 


I 4po _ 

3ml/ 


-2 


4 


(14) 


where the last approximation applies if po/pit) ^ 1 where po is the energy density at 
some early time origin. This gives us the relation between t and T 


3 ITOmpi 

271 V p* 


0.14 


mpi 

rp2 


(15) 


and At = tc{T^/Tl — 1) ~ 0.4 tc- Using ([l5|) evaluated at and solving (|T^ numerically 
with K = 1/2, ttw = 0.0336, B = 1.94, = 130 GeV, we get the sphaleron bound: 


C(Tc)>45.5 


c/iTc) 

Tc 


> 1-21 


(16) 
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where the second inequality was obtained using the fact that C(^c) is related to the order 
parameter (j){Tc) of the phase transition by: 


C(TJ = 

V c 

We can now check that with a more accurate evaluation of the integral in 
(0) is modihed intoQ 

In (1^) + 6 In era - era + Mra) seo 

weakening the sphaleron bound slightly: ^42.9 —> 4>c/Tc ^ 1.14. 


(17) 

), equation 

(18) 


2.1 Particle physics constraints from the sphaleron bound. 


As stated in the introduction, the relation ([TB|) is usually translated into a bound on the 
higgs mass. Indeed, 0(Tc)/Tc can be computed precisely from an analysis of the effective 
potential of the higgs at high temperature. For instance, a simple perturbative evaluation 
at one loop gives (with no daisy resummation) in the Standard Model 



SM 


2(2m^ + m|) 
3m| X nv 


(19) 


Plugging this expression into (p!6D leads to rrih < 37 GeV. As mentionned in the introduc¬ 
tion, more precise calculations show that in the standard Model [Q there is no hrst order 
phase transition at all for higgs masses rrih > 72 GeV. To make the phase transition hrst 
order, new bosonic degrees of freedom with large couplings to the higgs are needed. The 
MSSM precisely provides these, in particular, the right stop Ir with large Yukawa coupling 
can signihcantly enhance the strength of the phase transition [0,^,^ . Next formula 
(pop is an expression for 0(Tc)/Tc; it is not rigorously correct because it corresponds to a 
one-loop calculation and also to the ideal special case where 0 which means 

that the high temperature correction to the right stop mass is compensated by its 
negative soft mass parameter However, it is helpful to get a havor of what controls 
the strength of the EWPT in the MSSM: 




3/2 


\ T 
\ ^ / 


MSSM 


V Tc J 


nvm\ 


( 20 ) 


if |At|<mQ 


2 / = P VB{t')dt' = 2M.r] f°°dCC^e ^ where we treated MT) as a constant iMT) changes 

dowly with T), such approximation overestimates slightly the wash out. We noted G = C,{Tc)- With 
dCC'e-^ = ELo (G > 10), we get / = = 2VB{ta)tclCc. 
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Even if this one-loop evalnation nnderestimates (f){Tc)/Tc. compared with lattice calculation 
we see from this expression that the condition of preservation of the baryon asymmetry 
sets serious constraints on rrih, At/mq {mq is the soft mass parameter of the left stop 
ti) as well as mfj (and therefore ^ 1 ^). It is also instructive to keep in mind the one-loop 
expression for the Higgs mass (in the usual limit where the mass rriA of the CP-odd neutral 
Higgs boson is large): 


m, 


= cos^ 


2(3 + 


ml 




log 


9 2 

mA mA 

tR I L 

mf 


1 + 0 



( 21 ) 


While small values of At/mq < 250/1000 are required for baryogenesis, large values of 
At/mq > 250/1000 are required to satisfy the experimental bound on the higgs massQ 
and are in conflict with the sphaleron bound. 

The only alternative to reopen the window for electroweak baryogenesis is to weaken 
the sphaleron bound by modifying the last term in the left hand side of inequality (p^Sf ). 
For instance, to allow higgs masses in the range 115 GeV < mh ^118 GeV with large 
At/mq 700/1000 and unconstrained right stop mass (he. positive mfj values) one 
typically gets 0(Tc)/Tc 0.3 corresponding to C(^c) ~ 11.3. Therefore one would need 
to weaken the sphaleron bound by a factor 4. 


2.2 Evading the sphaleron bound. 


Let us find out what value for the expansion rate would be needed at a temperature 
Te 130 GeV to preserve the baryon asymmetry and be consistent with the larger window 
for the higgs and stop masses as presented below he. when using the value C(^c) ~ 11- 
Replacing C{Tc) ~ 11 into dill) leads to 


In(TcAf) < 10.1 


( 22 ) 


In a radiation dominated era with the relation (pAf) we get 


0.4T, 


H >. 10“’GeV 


(23) 


while we have if ~ 3 10 GeV so H is too small by 11 orders of magnitude at the 
EWPT. 

Earlier proposals to weaken the sphaleron bound have been to relax the standard as¬ 
sumption of radiation domination at the EWPT and assume instead a universe dominated 
by the energy density of a scalar field p0|-^. The main observational constraint is that 
such domination must terminate by the nucleosynthesis epoch. In p0|,^, it was shown 


that the energy density in a kinetic mode of a scalar field (which scales like a{t) faster 
than in the radiation case a(f)“"^) can significantly weaken the sphaleron bound and that 

^An ultimate alternative to both satisfy the bound on the higgs mass and be consistent with small 
values of AtjmQ would be to push the left stop mass to 2 or 3 TeV instead of 1 TeV. 
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such a phase of “kinetion” is required to evade large entropy release diluting the asym¬ 
metry. In it was argued that a universe dominated by the energy density of a more 
conventional scalar field (not necessarily in a kinetic mode) could also evade the erasure 
condition if the reheat temperature was as low as ~ MeV. However, the baryogenesis 
mechanism would have to be much more efficient than in the standard case thus standard 
electroweak baryogenesis does not work in models with low reheat temperature and one 
would need to implement an Affleck-Dine mechanism for instance. 

In this paper, we do not get rid of the radiation domination assumption and therefore 
do not have to assume any particular behaviour of a scalar field driving the expansion 
of the universe before nucleosynthesis. Our point of view is quite different. We are 
instead allowing deviations to the Friedmann equation itself at early times and study 
its consequences for electroweak baryogenesis. Let us assume for example that before 
nucleosynthesis the expansion rate behaved like 


H ~ 


piT) 

M3 


(24) 


where M is some scale to be discussed later. Such behaviour may a priori seem ad hoc 
but will be motivated in the next section. For a radiation dominated universe we get 


H ~ 


30 M3 


(25) 


leading to 


pit) ~ po 



4 Po t 
M3 


-1 


and the relation between time and temperature is now 


t ~ 


y \g*iTo)To 


\ 15 m3 _ 15 m3 

j 2^1^ 2^1^ 


(26) 


(27) 


so that At = tc — 1^ ~ 0.95 tc- Equation ( p2[) now translates into M < 20 TeV. More¬ 
over, for M < 15 TeV, the condition (|^ for the preservation of the baryon asymmetry 
is always satisfied whatever the value of CcQ. Before concluding that the expansion rate 
(^) can successfully lower the sphaleron bound or even lead to the disappearance of that 
bound for sufficiently low scale M we have to check two more things: 1) that sphalerons 
are at equilibrium in the symmetric phase i.e. that baryon number violation is still com- 
pletly efficient outside the bubble, 2) that electroweak interactions are at equilibrium as 
well so that usual techniques to compute the baryon asymmetry continue to apply even in 
such rapidly expanding universe. 2) is actually automatically satisfied when 1) is satisfied 
since the sphaleron processes are the slowliest ones, being suppressed by a high power of 

"^Weaker bounds on M are obtained when evaluating more properly the integral ®) in the same way 
as equation was obtained. 









aw] so we have to compare the rate of baryon number violation in the symmetric phase 
with the expansion rate H and require: 

M > H (28) 

According to latest calculations 0: 

rf ™ = (25.4 ± 2)alT^ (29) 


Using the standard expansion rate of the universe one obtains from (p8D the condition 
T < 7.10^^ GeV for sphaleron to be at equilibrium in the symmetric phase. On the 
other hand, for M ~ 20 TeV using the expansion rate (P^ ) we get T < 180 GeV. This 
temperature is now close but still above the critical temperature. It is actually sensitive 
to the exact prefactor in (p^) as well as the value of the numerical factor in the right 
hand side of equation (|2^) (which we took 0{1) in a first attempt). So a more accurate 
analysis will be needed to determine which precise values of M and the prefactor in H 
are consistent with sphalerons at equilibrium in the symmetric phase. 

One might also worry that in a fastly expanding universe it becomes difficult to nu¬ 
cleate bubbles and that the nucleation temperature T„ might be quite different from the 
critical temperature Tc at which the effective higgs potential has two degenerate minima. 
However, it has been shown that the bubble nucleation temperature depends only 
weakly on the expansion rate of the universe so that ~ T„ remains valid even when 
varying the expansion rate H over orders of magnitude. Gonsequently, the velocity of the 
phase interface may be subject to minor changes only. The same remark applies for other 
quantities which determine the baryon asymmetry generated like the wall thickness and 
perturbations in thermal population densities of particles at the vicinity of the wall. 

From this rough analysis we remark that one can substantially weaken the sphaleron 
bound using the expansion rate (^) as far as the scale M is low enough. Let us now 
study in which context this kind of behaviour might arise and how much freedom we have 
with respect to observational constraints. 


3 Two motivations to study deviations from Fried¬ 
mann equation 

1) No probe of FRW cosmology before nucleosynthesis 


Friedmann’s equation relates the expansion rate H of the universe with its energy 
density p{T) at a given temperature T. It is derived from the 00 component of Einstein’s 
equations in four dimensions, applied to the Robertson-Walker metric 

=—dt‘^ + a^{t){- - j^+r'^dQ'^) (30) 
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where a{t) is the scale factor and k the curvature parameter^ and reads: 


= 


a\^ p{T) k 


aj 


( 31 ) 


‘'PI 


Equation (|^ is a key ingredient of the standard model of cosmology, it serves as a basis 
to compute freeze-out temperatures and the density of relic particles in the universe. 
One of its major test comes from the success of the theory of primordial nucleosynthesis. 
However, there is no observational test of the Friedmann equation at times earlier than 
the primordial nucleosynthesis epoch characterized by an energy density p ~ 0.1 MeV^, 
and it is legitimate to question the extrapolation of FRW cosmological evolution before 
that epoch. Therefore, one could, from an empirical point of view, allow deviations to 
the expansion rate of the universe of the form: 


= 


p{T) , . .. 


(32) 


PI 


where cx{T) is a monotonic increasing function of temperature with a constrained value 


at nucleosynthesis, aiT = Tnud) ^1/10 


The origin of these deviations may come 
from modihcations to the Robertson-Walker ansatz (PDD by for example questioning the 
homogeneity and isotropy assumptions or it may come from allowing modihcations to 
general relativity itself as suggested by string theory. The existence of extra dimensions 
may induce corrections to the standard four-dimensional gravity, which are undetectable 
at low energy but might have played a signihcant role in the early universe. 


2) Deviations to Friedmann equation arise in higher dimensional theories of gravity 


The parametrization (p^ is actually motivated by results obtained in the context of 
brane cosmology. Any higher dimensional theory of gravity which reproduces 4d Einstein 
gravity is expected to reproduce Friedmann-Robertson-Walker cosmology at least at late 
times. One famous example is the 5d Randall-Sundrum (RS) model where the warped 

non factorizable geometry allows 4d Einstein gravity to be recovered on the 3-brane. For 
instance, the effective non-relativistic gravitation force between two masses mi and m 2 on 
the brane obeys the usual Newton’s potential plus some corrections scaling like 0{r~^): 


V{r) = 14ewton(r) (1 -h 


kPr‘^ 




.W = 


mim2 1 


m 


(33) 


PI 


where k is the inverse of the AdS radius of the RS geometry. The validity of the Newton’s 
formula has been checked down to distances r < 200 pm so that sets the constraint 


k >, IQ-^GeV 


(34) 


The cosmology of the RS model has been studied by perturbing the static solution of the 
5d Einstein equation and adding some energy density on the brane responsible for the 


^In the rest of the paper we will set k = 0. 
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expansion p^,|38|. It was shown that the expansion rate of the nniverse is of the form 
(^) with 


a{T) = O 


p{T)m 


PI 


M6 


(35) 


where M is the hve dimensional Planck mass. It is remarkable that such a modihed 
Friedmann equation is not in conflict with primordial nucleosynthesis provided that M > 
10 TeV while it leads to drastic changes at early times for which a{T) dominates and 
therefore cx instead of the conventional behaviour cx p |]3^,|40[|. It turns out that a 
more severe constraint on M comes from (^) which can be translated into M > 4 10"^ TeV 


by using the relation between the 5d and 4d Planck mass in the RS model = kMpi. 
Such bound on M is in conflict with our aim to weaken the sphaleron bound^ However, 
the result (|32D with (|35| ) can be generalized to any codimension one brane universe (for 
instance with vanishing tension and vanishing cosmological constant in the bulk i. e. for a 
relatively small warping) with compact stabilized hfth dimension. It has been extensively 
discussed in that provided the 55 component of Einstein equation has been used to 
stabilize the radion one generically recovers FRW cosmology for any codimension one 
brane modelQ. Let us consider the particular case of an Arkani-Hamed-Dimopoulos- 
Dvali type of geometry [Q. As is well-known, in the case with one extra dimension 
the fundamental Planck scale M cannot be lowered to the TeV scale. This would lead to 
major modihcations of gravity at observable distances. Therefore, if we want M ~ 10 —50 
TeV as required to lower the sphaleron bound we need at least more than one large extra 
dimension. So our point is to extrapolate the behaviour (|3^) to brane models with a 


higher number n of extra dimensions where M can be lowered to a scale of order O{10 
TeV). We do not know how the Friedmann equation looks like in the general case with 
n > 1. While brane cosmology has been the subject of intense investigation in the 
last two years it is still at its infancy. Most studies have focused on the 5d case and 
also were essentially concerned with the recovery of the FRW cosmology at low energy. 
An important point to stress is that the result (|3^j35|) is obtained hrst once the radion 
is stabilized and second in a perturbative approach by linearizing Einstein equations 
with respect to the perturbation 6a(t,y) (in a(t,y) = a{y)Rs + ^ci(t,y)) and p. What 
we are interested in in this paper is not the recovery of FRW but precisely the regime 
where this perturbative expansion breaks down and the nature of the modihcation to 
the usual Friedmann equation. A specihc feature of brane cosmology in 5d are the p^ 
term corrections. They are somehow a consequence of the junction conditions on the 
brane ||^ which are derived by identifying the singularity of the second derivative of 
the metric a”/a with the singularity of the energy-momentum tensor. To determine the 

® Also, if one wants to use the RS geometry to solve the hierarchy problem then two branes are needed 
and M mpi which is much too large to evade the sphaleron bound. 

^In brane-world geometry one can find solutions to Einstein’s equations where our four common 
dimensions only are expanding while the fifth extra dimension can be stabilized. This is to be contrasted 
with Kaluza-Klein cosmology where the energy density is uniformly distributed in all dimensions and leads 
to the expansion of the compact fifth dimension as well, consequently to a problematic time-dependant 
Newton ‘constant’. 
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junction conditions in the higher dimensional case one would need to solve equations 
of the type Aa(?/i, /a p6"‘{yi, ...,yn) where A is the laplacian in n dimensions. 

Being ignorant about the non static solutions of Einstein equations in the n > 1 case we 
will assume that terms still arise in Friedmann equation and are suppressed by as 
in the 5d case so that they dominate over the usual term at early times and in particular 
at the time of the EWPT. We will therefore use the expression 


H = a 


P 

M3 


(36) 


where a parametrizes our ignorance on the exact numerical factor; we will focus on a 
typical range of values [0.1 — 1]. 


4 Some comments on the cosmology with non stan¬ 
dard Friedmann equation 

Note that the transition between the standard and non standard cosmologies would take 
place when ~ ap/M^. In a radiation dominated era this leads to 


trans I - - I 

\mpi 


(37) 


With a = 1/5, M = 20 TeV we get ~ 1.4 MeV i.e. p ^ 135 MeV^ ~ 1350 
Therefore, transition to conventional cosmology takes place slightly before nucleosynthe¬ 
sis. 

While it is reasonable to expect that the highest temperature in the universe T,.eheat 
was 0{M) it has been argued that in theories with low fundamental Planck scale T,.eheat 
has to be much lower (< 1 GeV) to be compatible with observational constraints P5|,^ 


This is due to the presence of light modes other than standard model particles which can 
propagate in extra dimensions, in particular the graviton, and can affect substantially 
the energy density of the radiation after the decay of the inflaton on our brane. There 
are different independant kinds of bounds on T^eheat- They come from requiring that 1) 
the energy density of Kaluza-Klein gravitons does not modify the expansion rate of the 
universe at the time of nucleosynthesis, 2) the decay of KK gravitons into photons does 
not generate distorsions in the 2.1 K photon background radiation, 3) the energy density 
of KK gravitons does not lead to overclosure of the universe, 4) the cooling of the energy 
density of the radiation on our brane is dominated by the expansion and not the escape of 
gravitons into the bulk. All these constraints are based on the assumption that the bulk 
is empty (and static) and gravitons can only decay on our brane. Therefore they can be 
evaded if assuming that graviton can decay on other brane than ours. Also, a toroidal 
geometry is assumed for extra dimensions. Other choices of geometry lead to different 
spectra of KK modes which lead to weaker bounds on T,.eheat [0 and the occurrence of 
the EWPT in a radiation dominated era is no more questionned^. 

®For other discussions on the reheating temperature in non conventional brane cosmologies see also 


12 









It is interesting however to come back on condition 4). In our previous analysis of 
subsection 2.2, the relation between time and temperature was derived assuming that 
cooling is due to the expansion of the universe so that p{t) ~ and p|exp = —^Hp. 

The cooling due to the escape of gravitons into the extra dimensions was evaluated in ||43||: 


pU ~ -r”+VA/”+" 


(38) 


Requiring that the cooling by escape of graviton is less than ten percents leads to Hp > 
/2. For a standard cosmology this gives 


rjnn-\-l 


2 


30 ) mpi 


(39) 


If M = 10 TeV we obtain T < 13 GeV for n = 2. On the other hand using the expansion 
rate (^) we get 

T”-'s, (((^) M"-' (40) 

So, even if we were to assume the validity of (l38|) , at any temperature below M, the 
cooling is well dominated by the expansion. 

Let us hnish this section with the evaluation on the lower bound on M obtained by 
requiring that the expansion rate (pO]) should be suppressed at the time of nucleosynthesis 
when p ~ 0.1 MeV^. Requiring a{T) = < 1/10 leads to 


16.1 TeV 


(41) 


so that for cr = 1, M > 16.1 TeV, for a = 1/5, M > 9.4 TeV and for a = 1/10, M > 7.4 
TeV. 

In the next section we rehne the rough analysis of subsection 2.2. Our aim is to study 
the parameter space leading to a significant weakening of the sphaleron bound. There 
are three parameters to work with; The higher dimensional Planck mass M, the factor a 
appearing in the expansion rate (^) and k the prefactor in the sphaleron rate (^. We 
will limit our study to values of k in the range [0.1 — 1] since non perturbative calculations 
of the sphaleron rate indicate n >0.1. 


5 Condition to preserve the baryon asymmetry using 
non conventional cosmology 

Let us now analyse in more details the condition S > 10“^ expressing the preserva¬ 
tion of the baryon asymmetry produced at the EWPT. The integral in (^ is now I k. 
4:M.T]TctcCc^~‘’^, leading to 

In (j)))/]) + 61nC(ro) - m + HTcQ SO (42) 
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which is very similar to 


It can be rewritten as 




(43) 


where 


f{^,M,K,a) = 3.06+ ln( 

C 


K 


M\ 


aT^ 


+ 61n(^)- 37.52 ^ 

C ^ C 


(44) 


This last expression has been obtained by replacing in (|^ by its expression (|27|) and 
that is what makes the difference with the conventional cosmology. In figure |I| we plotted, 
for different values of M, /(|r, M, k, a) as a function of |r in the case a = 1, k = 0.2. 
We find that for M <23.1 TeV, / is always negative which means that the condition of 
preservation of the baryon asymmetry is always satisfied whatever the value for cj)c/Tc. 
Therefore, if a = 1 and M is in the range [16.1 — 23.1] TeV there is no sphaleron bound. 
For larger values of M, / becomes positive for certain values of ipc/Tc and therefore there 
is a sphaleron bound. Notice that M = 23.1 TeV corresponds to a discontinuity in the 
bound on (pc/Tc- From no constraint on (pc/Tc we jump to the constraint pc/Tc > 0.16. 

In figure 0, we summarize the status of the sphaleron bound in different regions of 
the parameter space {M,k) for the two cases a = 1 and a = 1/10. The horizontal 
line corresponds to the lower bound on M dictated by compatibility with nucleosynthesis. 
Another type of lower bound on M is obtained by requiring that sphalerons are at thermal 
equilibrium in the symmetric phase (see equation (^)) at Tc- It leads to the condition 
M > 0.1 Tc For = 130 GeV we get M > 14.3 a^/^TeV, which is slightly weaker 
than the bound obtained from nucleosynthesis. Anyhow, it is interesting to point out 
that independently from the nucleosynthesis constraint, what electroweak baryogenesis 
singles out is a window of a few tens of TeV for M. The upper bound corresponds to an 
expansion rate which is too small to freeze the baryon asymmetry in the broken phase 
at the EWPT while the lower bound corresponds to an expansion rate which is so large 
that sphalerons become out of equilibrium in the symmetric phase as well. The upper 
curves corresponds to pc/Tc = 0.3 which we presented in section 2 as the value obtained 
for higgs masses as large as 115 GeV and stop masses no more constrained to be lighter 
than the top mass. The lower curve delimits the regions where there is a sphaleron bound 
at least pc/T^ >0.15, and no sphaleron bound at all. 

We should add a last clarifying remark here. The main formula we have made use 
of in our analysis is equation (^) which gives the rate of baryon number violation in the 
broken phase. One should keep in mind that this formula is valid in a relatively narrow 
temperature interval Mw(T) « T « Esp(T). As stated in ||26|, for Egp/T < 1 
i.e. for p/T < 0.1, configurations other than the sphaleron contribute to the probability 
of making a transition with changing baryon number. Gonsequently, we can no longer 


apply the analysis of decay of metastable states with allowance for only the saddle 
point {i.e. sphaleron solution) of the energy functional. As p/T reaches zero, the sphaleron 
rate is expected to approach its value in the symmetric phase and the turnover around 
P/T 0.15 in the graph of figure 1 looses meaning. This however does not alter the 
validity of our qualitative results. Moreover, in the range p/T ~ 0.2 — 0.4, formula (^) 
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leads to an overestimation of the valne of the sphaleron rate (it slightly overpasses the 
valne in the symmetric phase) which means that by nsing that formnla in that regime we 
have slightly nnderestimated the freeze ont. 


6 Conclusion 


In snmmary, we stndied the condition for preservation of the baryon asymmetry prodnced 
at the EWPT when qnestioning the nsnal assnmption on the thermal history of the 
nniverse. To weaken the sphaleron bonnd by a factor 3 or 4 as reqnired to allow higgs 
masses as large as 115 GeV as well as a wider window of the MSSM parameter space, 
the expansion rate of the nniverse at the EWPT has to be increased dramatically (by ten 
orders of magnitnde at least) in order to freeze ont the sphaleron processes in the SU (2) 
broken phase. On the other hand, the expansion rate shonld not be larger than a critical 
valne above which sphalerons start to be inoperative in the symmetric phase as well. It 


tnrns ont that the the expansion rate (|36|) inspired by brane cosmology possesses these 
properties for M in the approximate range 15 — 50 TeV depending on the precise valnes 
of the parameters k and a. It wonld be interesting to check that the Friedmann eqnation 
we have postnlated effectively arise in generic brane models and more generally to stndy 
deviations to Friedmann eqnation before the nncleosynthesis epoch in the spirit of 
which do not necessarily rely on the existence of extra dimensions. 

The possibility we have investigated is one among others. For instance, we are still 
assnming radiation domination. One could imagine various ways to modify the cosmology 
at the EWPTQ. Also, in our work, M is the fundamental Planck scale, but it could corre¬ 
spond to a different scale. For example, the authors of found that the corrections 
in Friedmann equation could be suppressed by powers of Kw (in a matter dominated era) 
where Kw is a scale associated with the mass of the radion. Note that we are considering 
a regime where gravity becomes essentially higher dimensional and where radion modes 
could play a role. One thing to care about is that decay of these modes will not dilute 
the baryon asymmetry. 

Finally, we have focused here on the “out of equilibrium” condition of electroweak 
baryogenesis. There is still much activity in the electroweak baryogenesis community to 
determine whether there are enough CP violation sources in the MSSM to account for 
the amount of baryon asymmetry observed today {hb/ s ~ 10“^^). 
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f((^/T,M,K=0.2,a=l) 



Figure 1: /(|r, M, k, a) expresses the condition of preservation of the baryon asymmetry. 
For negative values of /, baryon asymmetry is preserved. This leads to a bound on the 
value of (pc/Tc for M > 23 TeV which depends on the fundamental Planck scale. As stated 
at the end of section 5, one should not pay attention to what happens to the graph for 
4^c/Tc <0.15. The turnover is an artifact of the breakdown of the approximation leading 
to formula (H). 
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Figure 2: Sphaleron bound in the parameter space {k,M). In Region I, there is no 
condition at all on (j)c/Tc- In region II, 0.16 < (j)c/Tc < 0.3. In region III, (j)c/Tc > 0.3. 
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